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ABSTRACT 

We describe improved modelling of the emission by dust in a toroidal-like 
structure heated by a central illuminating source within Active Galactic Nu- 
clei (AGN). We chose a simple but realistic torus geometry, a flared disc, and 
a dust grain distribution function including a full range of grain sizes. The 
optical depth within the torus is computed in detail taking into account the 
different sublimation temperatures of the silicate and graphite grains, which 
solves previously reported inconsistencies in the silicate emission feature in 
type-1 AGN. We exploit this model to study the spectral energy distributions 
(SEDs) of 58 extragalactic (both type-1 and type-2) sources using archival 
optical and infrared (IR) data. We find that both AGN and starburst con- 
tributions are often required to reproduce the observed SEDs, although in 
a few cases they are very well fitted by a pure AGN component. The AGN 
contribution to the far-IR luminosity is found to be higher in type-1 sources, 
with all the type-2 requiring a substantial contribution from a circum-nuclear 
starburst. Our results appear in agreement with the AGN Unified Scheme, 
since the distributions of key parameters of the torus models turn out to be 
compatible for type-1 and type-2 AGN. Further support to the unification 
concept comes from comparison with medium-resolution IR spectra of type- 
1 AGN by the Spitzer observatory, showing evidence for a moderate silicate 
emission around 10 pm which our code reproduces. From our analysis we in- 
fer accretion flows in the inner nucleus of local AGN characterized by high 
equatorial optical depths {Ay ~ 100), moderate sizes {Rmax < 100 pc) and 
very high covering factors (/ ~ 80 per cent) on average. 
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nus - galaxies: individual: NGC 1068 - galaxies: individual: Mrk 231 - galax- 
ies: starburst - infrared: general 



1 INTRODUCTION 

The idea that the quasar infrared (IR) spectral energy distribut ion (SEP) could be dorn - 



Neugebauer et al. 



(19791), 



inated by thermal emission from dust was already proposed by 
who noted an excess at ~ 3.5 pm with respect to the po wer^law emission that was known 
to dominate the ultraviolet (UV) and optical spectrum. iBarvainisI ()l987l ) was the first to 
demonstrate that it was possible to describe this emission by means of thermal radiation 
from dust heated by the primary optical/UV continuum known to dominate type-1 AGN 
SED. He considered an optically thick, either smooth or a clumpy dust distribution within a 
non-spherical geometry. Such a model could naturally account for the observed 3 pm bump, 
which was attributed to graphite grains with temperatures close to their sublimation limit. 
The broadness of the IR SED, that could not be reproduced by means of a single-temperature 
black body, was explained in terms of multiple temperature components of dust. 

More quantitative models have been developed since then, accounting for different fea- 
tures and using diff e rent c omputational approaches for solving the radiative transfer equa- 



tion. iPier fc KrolikI (|l992i ) suggested an emission model from an annular ring of dust. As- 
suming a common extinction cu rve for all of t he gr ains, a distribution function of grain 



composition and sizes as given in 



Mathis et al, 



()l977f l and not accounting for light scatter- 



ing, they solved the radiative transfer problem with the Bohm-Vitense method described in 
Mihalad (|l97g l. Based on this model, they explored the main properties of IR emission in 
AGN, such as the width of the SED and the wavelength of peak emission, IR colours and 
the presence and /or absence of the 9.7 pm silicate feature. 



Granato fc Danesd ()l994[ ) solved the radiative transfe r equation bv me a ns of the A 



iteration method (see 12.61 for more details), adopting the iRowan-RobinsonI ()l992n model 
for interstellar dust, based on both graphite and silicate grains, and taking into account 
radiation scattering. Their adopted geometry was the so-called flared disc, with dust grain 



de nsity that w'as al 



owed to vary both along the radial and the vertical coordinates. 



Stenholml ()l994r ) studied a configuration given by an optically thick silicate disc, finding 



no need for silicate grains depletion or non-standard dust composition in order to weaken 
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the strength of the 9.7 pm feature, which turned out to be weak in most of his models 



Pier&Krolik 



with o ptical depth larger than 1. Furthermore he found that, in contrast with 
( 19921 ) ■ scattering is not negligible and does affect the shape of the near IR emission in those 
mo dels where optical depth is high . 



Efstathiou fc Rowan-RobinsonI ()l995[ ) analysed three different geometrical configurations 



for the torus, a flared and a tapered disc and an anisotropic sphere, and found that a tapered 
disc (i.e. a disc whose height increases along with distance from the center but flattens at 
a constant height in the outer parts) with an opening angle of 45° was the most succesful 
representation of what was o bserved. They followed a. ray-tracing numerical method, with 
the same dust mixture as in iRowan-RobinsonI ()l992| ). The sublimation temperature was 



assumed to be 1000 K for all grains, while the distances from the central source where 
this temperature is reach ed were assumed to be dependent on the grain size. In the same 



line, 



Manske et al 



( 1998f ) modelled the dusty torus as a flared disc in which both graphite 
and silicate where taken to sublimate at a temperature of 1500 K. In their analysis they 
concluded that flared tori, dust composed of silicate and graphite grains and a power-law 
density distrib ution failed to explain the observed featureless AGN spectra. 



In parallel, 



Nenkova et al. 



(2002) developed a torus model in which dust was distributed 



in clum ps. Thev solved the r adiative transfer for each clump by means of the radiative code 



DUSTY ( Nenkova et al 



19991 ) and then computed the final emission as the sum of clumps 
emission calculated using a ray-tracing approach. 



At variance with 



Efstathiou fc Rowan- RobinsonI (|1995| ) . Ivan Bemmel fc Dullemondl (|2003[ ) 



claim that the obscuring torus could not be flat (tapered, for example), and should rather 
be conical or flaring. Their more reliable models predicted constant density over the whole 
torus and a quite large value for the inner radius (~ 10 pc), which they found necessary in 
order to fit observed colours. They also trie d to modifv th e powe r-law distribution function 



for grain sizes with respect to that given by 



Mathis et al, 



( 19771 ). in order to attenuate the 



strength of the 10 pm sili cate feature and recently conclud ed that there was no need for a 



clumpy dust distribution (|Dullemond fc van Bemmelll2005[ ) 



In the present paper we try to combine most of the successful features of the previous 
works to built an improved model explaining IR emission in AGN. The paper is organised 
as follows. Section El presents the physical and computational details of the model. Sectional 
discusses the contribution of each of the model parameters in the produced SEDs. In Section 
lUwe explain the observed SEDs of a series of objects taken from the literature using various 
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Figure 1. Shape of the toroidal dust distribution in our model. In this figure R is the ratio of the external to internal radius 
and © is the full— opening angle of the torus. 

realisations of the torus model in combination with starburst templates. Finally in Section 
El we discuss our results and conclusions. 



2 MODEL DESCRIPTION 



2.1 Torus Geometry and Grain Sublimation 

We use a simple but realistic geomet r y for the spatial dist ribut ion of dust: the flared disc 



1995 



Manske et al. 



199S and 



van Bemmel fc DuUemond 



f Efstathiou fc Rowan- Robinson 

rl 

120031 ). that can be represented as two concentric spheres, delimiting respectively the in- 
ner and the outer torus radius, having the polar cones removed. 

The size of the torus is defined by the outer radius Rmax ^ the inner radius being defined 
by the sublimation temperature of dust grains under the influence of the strong nuclear 
radiation field - and by the angular opening angle 6 of the torus itself (see Fig. . Consid- 
ering a typical grap hite grain with ra dius = 0.05 pm and using the thermal equilibrium 
equation (see Eq. I7jl . iBarvainisI ()l987l ) derives the following formula to compute the distance 
at which this grain reaches the sublimation temperature: 

Rm^n^l■^■^|L^ -T^^,^ M, (l) 

where Lf^^ is the bolometric ultraviolet/optical luminosity emitted by the central source, 
expressed in units of 10^^ erg s~^ and T1500 is the sublimation temperature of the dust grain 
given in units of 1500 K. Note that this relations depends on both the grain dimension and 
species but we here assume as a common minimum radius the average one. We account for 
the different sublimation temperature of silicate grains, however, yielding a minimum radius 
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Figure 2. The geometrical grid is defined with respect to a spherical coordinates system. The equatorial plane of the torus 
lies in the xy plane, so that when the view angle ^ = 90° the geometry reproduces a type-2 objects, while ^ = 0° would 
correspond to the type-1 view, according to the Unified Model. 

that can be up to three times larger than the one for graphite. Furthermore, the presence of 
the dust layer composed by graphite grains, that attenuates the X-ray and UV emission, is 
taken into account when computing the silicate minimum radius. 

2.2 Definition of the geometrical grid 

We define a geometrical grid with reference to a polar coordinates system (r, 6, (f), as shown 
in Fig. 121 by dividing the torus in volume elements. We then compute all the quantities of 
interest within a given volume element as if all dust properties were those of the center. 
Since the temperature gradient is higher in the inner parts, we used a logarithmic-spaced 
division for the radial coordinate, in order to better sample the properties of the emission 
in the innermost regions. 

The grid has a constant number of subdivisions in the azimuthal coordinate {(f) that 
is fixed to 36, while the equatorial distance coordinate is divided in 12, 14 and 16 parts 
for models with respectively 60°, 100° and 140° of torus amplitude. The radial coordinate 
is divided in 40, 80 and 120 parts for Rmax/Rmin = 30, 100 and 300 respectively. The 
wavelentgh values at which the spectrum is computed consist in a grid equally-spaced in 
logarithm with 120 points in the range from 10^'^ to 10^ jam. 

Neither the computation of the emitted radiation from each grain species in each volume 
elements within the torus, nor the computation of the optical depth between all torus cells 
are particularly CPU time consuming. However, for a typical ~ 15000 element grid (see 
Sect. 3.1 below) this computation must be repeated ~ 15000^ times, if we want to evaluate 
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the contribution to the radiation field in each element from all other in the torus. It is 
therefore essential to exploit all the symmetries of the adopted geometry, identifying groups 
of volume elements whose properties are identical all over the torus grid. For instance, 
all volume elements at any given ip have exactly the same dust density, temperature and 
emission conditions. It is sufficient to compute these quantities just for these elements, which 
we will refer to as sample elements. 

Furthermore, when computing the incoming radiation on a given sample element, the 
problem is symmetric with respect to the xz plane, so it is possible to compute incoming 
radiation from just one side of the torus. 



2.3 Dust properties 



The main dust components are silicate and graphite grains, in almost equal percentages. The 
former are responsible for the observed absorption feature at ~ 9.7 jam in type-2 objects, 
while the latter are responsible for the rapid decline of the emission at wavelength shortwards 
of a few Jim, corresponding to a black body emission of about 1500 K, the sublimation 
temperature for these particular grains. 

Most mod els in the literature make use of a small number of grairi categ ories at differ- 



ent sizes (e.g. 



Granato &: Danese 



approach similar to that of 



1994 : 



Efstathiou fc Rowan- Robinson 



van Bemmel fc Dullemond| " ()2003[ ). we decided to adopt a more 



1995^ . Following an 



complete distri bution using, fo r each of the two species, scattering and absorption coeffi- 



cients given bv iLaor fc Draind ()1993[ ) for different grains sizes from 0.005 to 0.25 pm and 
from 0.025 to 0.25 pm for graphite and silicate respectively. 

As for the distribut ion function for the grains sizes a we adopted the MRN distribution 



see 



Mathis et al 



19771) : 



dN{a) = 10^' • a'^da 



(2) 



where dN{a) is the number of grains with radius between a and a + da normalized to 
the number of hy drogen atoms. As f or the exponent q, the value of the standard galactic 
extinction curve (jMathis et al.l 119771 q = —3.5) has been adopted, so that the extinction 



curve that characterizes our model is similar to the galactic one. The constant Ai gives the 
normalization with respe ct to hydrogen abun dance, and it is taken to be —25.16 for graphite 



and —25.11 for silicates (Draine fc Lee 



19m . 
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2.4 The density law throughout the torus 

We adopted a law for the gas density within the torus allowing for a density gradient along 
the radial and the polar distance (6) coordinates: 

p(r,^) =«.r^.e-^^l^°^(')l. (3) 

The normalization constant a here is a linear function of the equatorial optical depth (see 
Section ITTj) . The latter will be used in the following instead of the a parameter. 



2.5 Spectral intensity of the central power source 

We assume that the torus is illuminated by a central point-like energy source with isotropic 
emission. We define its spectrum in the wavelength range 0.001 pm ^20 pm (corresponding 
to the frequency range ~ 3 • 10^''-^ ~ 1.5 • 10^'^ Hz) and described it as a c omposition of power 
l aws w ith variable indices . The values that we adopted are the same as in 



Jiggj) or 



Nenkova et al. 



Granato fc Danesf 



( 20021 ). namely: 



1.2 



A-L(A) = < 



Ln-A 



-0.5 



[erg/s] if 0.001 < A < 0.03 [fi] 
[erg/s] if 0.03 < A < 0.125 [fi] 
[erg/s] if 0.125 < A < 20.0 [fi] 



(4) 



where Lq sets the bolometric luminosity normalization constant. 



2.6 Solution of the radiative transfer equation 

The radiation emitted by the above described physical configuration is given by the sum of 
the primary source located in the torus center and a secondary contribution given by thermal 
and scattering dust emission. Therefore dust plays a double role, first absorbing the radiant 
energy emitted by the primary source, partially or totally obscuring it, then re-emitting it at 
longer wavelengths, typically in the spectral range 1 -r- 1000 \im. A complication arises from 
the fact that dust can be optically thick to its own radiation. This requires the solution of 
the radiative transfer equation [Eq. © below]. In particular the presence of a dust scattering 
component makes this equation not solvable analytically: the source function is not known 
a priori, depending on the solution Jj, itself. 

The thermal and scattered components of dust radiation are described by the absorption 
(oiy) and the scattering (cr^) coefficients respectively, so that the radiative transfer equation 
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will be written as: 



^ = - (a^ + cr^) ■ (/^ - S^) (5) 



in wh ich Si, is the source function, that can also be written as (see iRvbicki &: Lightman 



19791 ): 



oiiiBjj -\- (JiiJn , . 

Su = ■ , (6) 

a,, + a,y 

i.e. a weighted average of the two separate source functions. 

We solved the transfer equation numerically adopting the A-iteration method: A lower 
limit for the incident radiation is estimated assuming that the only source of radiation (and 
hence of dust-heating) is the central, non-thermal, source. The equilibrium temperature 
(and hence thermal emission) of the grains of each species in each sample element is found 
by solving the thermal equilibrium equation: 

/ QabsW ■ r\X)dX - / Qe^(A) • B{X, T,k)d\ = 0, (7) 

J Am J Am 

where the equation is written for a given species of grain and B{\,Tik) is the black body 
emission. Tj^ is the temperature of a given grain within the ik-th sample element, J"^^ repre- 
sents the incoming specific intensity on the volume element, written, for the first iteration, 
as follows: 



J'\\) = lAGN{tk,X) = — ■ -^exp[-r,fc(A)]. (8) 

An 47rr4 

Here rjfc(A) is the optical depth between the central source and the ik-th sample element 



and its distance. 

At the second iteration the new incoming radiant flux on each grain of a given geometrical 
element ik that will enter in the Eq. [7|of the thermal equilibrium will be: 



7V-1 

J^'W = IagnW + (9) 



e=l 
' Ndust 



J2 [47r.7ra?,Q7(A)fi(A,T,,,) + 



id 



where e indicates the generic volume element and is their total number. Jf^tW 
total specific intensity hitting the e-element and coming from both the AGN and the dust 
from the rest of the torus (scattering plus thermal) as computed from the previous step. 
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0.001 0.01 0.1 1 10 100 1000 0.001 0.01 0.1 1 10 100 1000 



Figure 3. Emission spectra by torus elements at increasing radial coordinate for dust torus distributions illuminated by a 
nuclear source of L^'^^ = lO**^ erg s~^. Two extreme cases are shown. On the left panel a model with equatorial optical depth 
r(9.7) = 0.1 and an outer to inner radii ratio of 30 is shown. The closest element to the central source is labelled 1, the most 
distant 40 (the different lines correspond to physical locations of 1.36, 1.91, 2.92, 4.50, 6.8, 10.4, 16.0, 24.4 and 37.4 pc). On 
the right panel the equatorial optical depth is taken to be t(9.7) = 10.0 and the outer-inner ratio is 300. 

These calculations are iteratively repeated until the difference in the temperature value 
between two steps falls below one Kelvin for each grain species in each volume element. If 
such small variations do not have big influence on the final result when dealing with relatively 
low-density environments, it is important to achieve a good convergence when dealing with 
high optical depth tori, to avoid losses of energy due to the numerical method adopted. 

The conservation of energy is always verified within 1 per cent for typical solutions and 
up to 10 per cent for solutions with very high optical depth and non-constant dust density. 
These discrepancies are due to some approximations in the maximum temperature of the 
grains, but do not affect any of our conclusions. 

In models with high dust densities where thermal emission strongly contributes to the 
self heating of the dust, the smaller grains can reach and even slightly exceed the sublimation 
temperature. In order to avoid further re-calculation of the optical depth and density profiles 
- that would result in a strong increasing of the computation time - temperatures which are 
too high are set equal to the sublimation value. If this also happens for the 0.05 pm grains, 
then the inner radius is moved further out with respect to the value computed by means of 
Eq.IH 

Various model solutions are reported in Fig. IHl where dust emission from elements at 
different distances from the central source are shown. The model solution on the left is built 
with a total opening angle of the torus of 140°, a ratio Rmax/ Rmin = 30 and a constant 
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density profile with equatorial optical depth req(9.7) = 0.1. The minimum radius turns out 
to be 1.30 pc for graphite and 6.27 pc for silicate grains. The contribution to emission from 
the silicate grains, which is easily recognizable from the feature in emission at 9.7 jim, only 
appears after the element number 15. The scattering component of radiation is clearly visible 
shortwards of A ~ 1 pm. 

The model shown on the right side of the panel has the same opening angle for the torus, 
an outer-to-inner ratio of 300, constant density profile and Tgg{9.7) = 10.0. In this case the 
minimum radius is the same (1.3 pc) for graphite, but it is lower for silicate grains (3.94 pc) 
whose emission appears after the element number 30 over a total of 120 radial subdivisions. 



2.7 Details on the computation of dust density and optical depth 

The model needs as a starting point the value of the equatorial optical depth at a given refer- 
ence wavelength. Starting from this input condition we calculate the gas and dust densities 
and the number of the different dust grains within each sample element. With reference 
to Eq. 121 providing the spatial variation of density through the (3 and 7 parameters, the 
problem is to find the a value yielding an equatorial optical depth Teg (Are/)- The expression 
for the latter is: 

rBmax (Ndust 

re,{Kef)= / <^ Yl [QtdiKef) + QU^ref)] (10) 

a-^TT ■ Ndid ■ pH{r,9)dr} 

where Ndid is the numeric density of the id-th grain, normalized to that of hydrogen. Taking 
this into account and the difference in the sublimation temperatures of graphite and silicate 
and the discontinuity in the chemical composition of dust along the radial coordinate, the 
density is a step function that can be written as: 

J . r < Rsii 

f^r\r,e)={ (11) 
[ pfr*(r,0) r>Rsu 

where the indices id refer only to the silicate component. An analogous condition applies to 
the graphite grain distribution and their corresponding id indices. Substituting the p from 
Eq. 121 in Eq. ^1 and solving for a we obtain the correct value of density and number of dust 
grains at each point of the grid. 

The most general expression for the optical depth between any two points -Pi(ri, 6i, ipi) 
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and P2{r2, O2, ^2), located within the torus is: 

/T2 PO2 ^^^^^ 
/ / E [^^'^(^) + Q'M X (12) 

7ra^^ Ndid ■ Pid'^r, 9) dip dO dr. 

The optical depth is calculated as the path-integral of the density function f{r,6,ip) = 
a ■ ■ e~'^^l™'^(^)l along the straight line between the two points. Considering the conditions 
im in IT^ automatically accounts for the fact that there are inner regions free of silicates due 
to their lower sublimation temperature and regions completely dust-free. 

2.8 Calculation of the emitted spectrum 

Adding the emission from all torus elements to that of the central source, attenuated by 
the correct value of optical depth when needed, yields the total emitted spectrum. Here it 
is computed for 10 different inclination values of line-of-sight, spanning from to 90° with 
respect to the equatorial plane. 

In Figs, m and El we show two emission models with different lines-of-sight and different 
physical and geometrical parameters. The first model is meant to represent a geometry in 
which the amplitude of the torus is 140°, the ratio between external and internal radii is 
Rmax/Rmin = 30, the density is constant and the optical depth is t(9.7) = 0.1. Such a 
configuration is obtained with a hydrogen column density of 9.0 x 10^^ cm^ and corresponds 
to an optical extinction of Ay = 2.3. Note that the shape of the SED remains quite constant 
with varying angle, as long as the central source is seen directly. A similar behaviour is found 
when both the optical depth and the torus size increase. The emission pattern of Fig. El has 
a radial ratio of 300 and the equatorial optical is t(9.7) = 10.0. 

3 MODEL REALISATIONS 

We have computed a grid of models by varying the different parameters. Some of their main 
features and a rough comparison with other models in literature are presented here. 

3.1 The grid of model solutions 

We obtained a set of models spanning the widest range of physical and geometrical char- 
acteristics of dusty tori discussed in the literature. The models have varying radial ratios 
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Figure 4. Emission spectra obtained for 10 different line-of-sight inclinations from 0° to 90° at regular steps of 10° . Here the 
torus amplitude angle is 140°, the ratio between the outer and inner radii is Rmax/Rmin = 30, the density is constant and the 
equatorial optical depth is t(9.7) = 0.1. 




Rmax/Rmin = 30, 100, 300, Corresponding to radial sizes from about 40 to about 400 pc 
(for an accretion luminosity of 10^^ erg s~^), while the torus full opening angles were taken 
to be 60°, 100° and 140°, corresponding to sizes of the torus between a minimum of ~ 20 
to a maximum of ~ 370 pc. The equatorial optical depth was always measured at A = 9.7 
pm and the models have been computed at values of t(9.7) = 0.1, 0.3, 0.6, 1.0, 2.0, 3.0, 
6.0 and 10.0, translating into values of absorption in the V band from few to few hundreds 
magnitudes, depending also on the Rmax/Rmin ratio. The parameters used to define the law 
for the spatial behaviour of density were j3 = 0, —0.5 and —1.0 and 7 = and 6 (see Eq. 
El for the definitions). In order to have a full exploration of the parameter space, we also 
considered models with density increasing towards the outer parts of the torus, obtained by 
setting (3 = 0.5. In Tabled we summarize the main features of the model grid. 
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Parameter 



Adopted values 




30, 100, 300 

0.1, 0.3, 0.6, 1.0, 2.0, 3.0, 6.0, 10.0 
0.5, 0.0, -0.5, -1.0 
0.0, 6.0 

60°, 100°, 140° 



7 

e 



Table 1. The parameters of the models of the grid. 



3.2 Main features of the models 

3.2.1 The distribution of grains temperatures 

One of the most evident features of the AGN IR emission is the presence of multiple tem- 
perature components. Especially when viewed face-on, the whole range of temperatures can 
be seen, from the innermost part - where grains are close to the sublimation limit - to the 
outer shells where the dust is colder, only heated by reprocessed IR radiation. 

For models in which the spatial profile of density is kept constant with respect to the 
altitude angle 6, the temperature is assumed to vary only along the radial coordinate. Models 
with non-constant density profile can show a temperature gradient also with respect to the 
altitude coordinate. In Fig. IHl we show the trend of the dust equilibrium temperature as a 
function of the coordinates of the sample element. Two extreme cases are shown: a model 
with t(9.7) = 0.1, Rmax/Rmin = 30 and a constant density, and a model with the same 
parameters but a density profile varying 

li]^Qr-^-0 .^-6x\cosie)\ (in both graphite grain 

with radius of 0.01 pm is considered). 

In the case of a constant dust density distribution, the temperature profiles drop roughly 
linearly with the distance from the central source. A non-constant density profile, in our case 
peaking in the inner part of the torus (close to the central source and to the equatorial plane), 
brings some differences. First of all, the decrease of dust density with increasing altitude on 
the equatorial plane results in lower values of the optical depth of the elements at lower 6 
coordinate. The temperature in the most distant volume elements will then increase with 
increasing altitude with respect to the equatorial plane. Furthermore, since the dust density 
peaks in the innermost regions, the temperature gradient will be higher for those volume 
elements that are closer to the inner edge of the torus. Moreover the high density (hence high 
optical depth) values that are found close to the inner radius are responsible for not always 
reaching the sublimation temperature: in such configurations the central source undergoes 
a high attenuation in the first layers of dust. 
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Figure 6. Temperature profiles of graphite grains with a radius of 0.01 iim as a function of the torus element. Each vertical 
blocks spans the radial distance of the torus and corresponds to a different polar angle. On the left panel a model with constant 
density profile, t = 0.1 and radii ratio of 30 is shown. In this case the only temperature gradient present is the one with respect 
to the radius. On the right panel a model with the same characteristics but with also a non-constant density profile (/3 = — 1.0 
and 7 = 6). The horizontal lines represent the values of sublimation temperature assumed for graphite (1500 K) and for silicate 
(1000 K) grains respectively. 



3.2.2 Dependence on the viewing angle 

As expected for a geometry with no spherical symmetry, the emission is highly anisotropic 
especially when the lines of sight that cross the torus and those that see the AGN directly 
are considered (see Figs. |3]and[S)). This feature is further highlighted if one considers a 
pattern in which the density is not constant along the altitude angle 6. 

Figs, m and El show that the SED remain quite constant even with varying line-of-sight 
angle, until the latter intersects the dust-free zone, i.e. at \1' < 90° — 0/2. The only part 
of the emission that changes is the one shortwards of 1 pm for low optical depth models, 
where the emission of the non-thermal central source is directly seen in the type-1 scheme 
configuration. In high optical depth systems this behaviour is maintained, but the drastic 
difference between type-2 and type-1 concerns the emission shortwards of 30 pm. 

If a non-constant density pattern is assumed the change is smoother since the emission 
from the central source becomes more visible as the torus is seen at decreasing polar angles, 
where the density is assumed to be lower with respect to the equatorial plane. Notice that 
with such density profiles and in those models with low optical depth values the direct 
view of the central source is achieved even when \E' > 90° — 0/2 since the optical depth at 
higher altitudes becomes small enough. Then a substantial contribution of the accretion disc 
emission can be seen even through the dust cone. Configurations with decreasing number of 
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dust grains with increasing altitude on the equatorial plane can be consid ered to mimic the 
tapered torus geometry analysed bv lEfstathiou fc Rowan-RobinsonI ()l995r ) since, at least for 
low values of the equatorial optical depth [i.e. t(9.7) < 1.0], the bulk of emission comes 
from those volume elements that are closer to the equatorial plane. The other elements of 
the geometrical grid will in fact provide a smaller contribution to the total emission, as they 
will contain a smaller number of grains. 



3.2.3 The width of IR spectrum 



One of the most important features of a dusty torus model is the capability of r eproducing the 



observed width of the IR spec t rum d ue to multi-temperature emission. Both 



and 



Pier fc Krolik 



Granato fc Danesd ()l994f ) use for its measure the logarithmic wavelength interval 



in which the power XF\ emitted in the IR is more than one third of the peak value (Wir). 
When measured on a black b ody this parameter has a v alue of ~ 0.7 while, for example, in 
the observed spectra used bv lGranato fc Danesd (jl99J) its value is always larger than 1.3. 
In Fig. [7|we report the trend of the bump width as a function of optical depth for different 
torus sizes, as viewed along the torus equatorial plane. This parameter depends both on 
the optical depth value and on the torus size: increasing the values of the optical depth will 
tend to generally shrink the bump since the inner, hotter, parts will become less visible. At 
rising values of the radii ratio the bump will tend to broaden, since larger radii imply higher 
amounts of colder dust. 

There are a few considerations to be made on the interpretation of the Wjr trend in Fig. 
[7| The maximum of emission for low values of optical depth falls right on top of the silicate 
feature which in these cases can displa y a very intense emiss ion. This will in turn give a low 
value for the bump width as defined bv lPier fc KrolikI (|l992| ). Furthermore, very high values 
of optical depth will cause the silicate feature to appear in strong absorption especially at 
equatorial optical depth higher than 6, resulting, again, in a low width for the IR bump. 

A density profile peaking towards the innermost regions of the torus is less affected by 
this behaviour. This is mainly due to the fact that in the latter case the silicate feature is 
less prominent both in emission and in absorption, compared to a constant dust distribution 
configuration. 
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Figure 7. The trend of the IR bump width is shown as a function of the optical depth and for various radii ratios. The 
parameter Win is defined as the logarithmic interval in wavelength, were the power emitted is more than one third of the 
maximum. On the left-side panel models with homegeneous density distributions are shown, while on the right panel a model 
with density parameters 7 = 6.0 and [3 = —0.5 is plotted. 



3.2.4 The peak of the IR emission 

Another important parameter that can be used to characterize a model solution is Xmax, the 
wavelength at which the IR emission - expressed as XF{X) - peaks. In Fig. |Hlthis parameter 
is shown against the optical depth and for the three torus sizes for which the model has 
been computed. Only type-2 (equatorial view) emission is shown here. The peak wavelength 
ranges from 9.7 pm, corresponding to the presence of a strong silicate emission in models 
with low optical depth, to about 60 pm obtained in the models with the largest outer radii 
and optical depth values. 

A non-constant dust density profile results in a hotter emission, with the emission peaking 
at smaller wavelengths. This occurs because lower densities - and hence lower optical depths 
- on the outer regions of the torus make the inner, hotter, regions more visible to the observer, 
increasing this way the average observed temperature. 



3.2.5 The 9.7 //m silicate feature 

The presence of the absorption feature at 9.7 pm in the spectra of Seyfert 2 galaxies has 
been used in the literature as a proof for the presence of silicate grains. At the same time, 
IR observations often showed rather featureless spectra for Seyfert 1 galaxies, which brought 
into question the reliability of the Unified Model: the characteristics of the silicate emission 



at this wavelength has become a crucial test of the validity of models. 



Nenkova et al. 



( 20021) 
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Figure 8. Dependence of the wavelength of the peak IR emission (expressed in XF^) on the optical depth and torus size for a 
constant (left panel) and a non-constant (right panel) density distribution in the torus, calculated for a type-2 viewing angle. 
The constant value of ~ 10 pm for small values of optical depth is due to the presence of the intense 9.7 pm silicate emission 
dominating the spectrum. Here the dependences on the optical depth and on the torus size are comparable. No significant 
difi'erences are found for the non-constant density model. 

suggested that a clumpy distribution of dust may provide an explanation for the lack of 



the silicate feature in type-1 sources, while 



Dullemond fc van Bemmell (|2005[ ) claim that the 



same result can be achieved in smooth tori. The depth of the same feature in type-2 sources, 
however, is shallower i n clumpy tori. 



Only very recently ISiebenmorgen et al.l (120051 ) and 



Hao et al, 



( 20051 ) have found clear 



observational evidences for the presence of the 9.7 pm silicate emission feature in type-1 
quasars from spectral observations by t he Spitzer Space Te l escop e. 



Following 



Pier fc KrolikI (|l992t) and 



Granato &: Danesd (|l994^ . we measured the ampli- 



tude of this feature by means of the dimensionless parameter 5*9,7, defined as the natural 

logarithm of the ratio between the actual flux emitted at 9.7 pm and the flux value at the 

same wavelength interpolated between 6.8 and 13.9 jim: 

F(9.7) 



S9.7 = In 



(13) 



LFc(9.7)_ 

In Fig.Elwe report the variation of the parameter 6*9.7 as a function of the external to internal 
radial ratios and of the optical depth for a type-2 (left panel) and type-1 lines-of-sight (right 
panel). A homogeneous distribution of dust is assumed for models in the upper panels while 
in the lower panels the behaviour of non-constant density models is shown. 

The strongest dependence of the silicate emission for type-2 viewing angles is on the 
optical depth, while the torus size plays a secondary role that becomes more important 
when a type-1 viewpoint is assumed. In the latter case the variations are much smoother 
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Figure 9. The variation of the 9.7 pm feature intensity 59.7 as a function of the optical depth for three torus sizes. In the 
two panels at the top, a torus with a constant density distribution is considered, while at the bottom a model with a density 
gradient is shown (/3 = —0.5 and 7 = 6.0). A type-2 line-of-sight = 90°) is considered on the left-side plots: here we see 
that the strongest dependence is on the optical depth rather than the torus dimensions. On the right-hand panels the type-1 
view ('5 = 0°) is shown. The variation of the S9.7 parameter is smoother in this case and the most extreme variations are 
encountered in models with the smaller sizes. Models with the density peaking in the innermost regions show a general smaller 
variation. 



and smaller torus sizes appear to be more sensitive to changes in optical depth. Hence, 
assuming a density profile that decreases towards the outer zones makes the interval of 5*9.7 
values more narrow with respect to a constant density profile, both in the type-2 and in the 
type-1 solutions. This is likely due to the fact that the silicate features in such models is 
more nested in high optical depth zones, and its presence is attenuated. 

Models in which the silicate feature is absent or even measured in absorption for a polar 
(type-1) view are obtained with high values of the equatorial optical depth [t(9.7) around 
10], smaller outer-to-inner radius ratio and 7 = 0. The amount of absorption increases as the 
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exponent /? decreases, while positive values yield an almost featureless spectrum. Models with 
decreasing density at increasing angular distance from the equatorial plane always display 
the silicate feature in emission, as the hotter innermost regions where the silicate feature is 
stronger are less extinguished when the torus is seen face on. The intensity of the Silicate 
emission feature predicted by our model - thanks to the detailed description of the inner 
torus boundary where the feature is mainly produced - is moderate, and almost absent for 
r(9.7) > 1 (see Fig. IHl bottom-right panel). 



3.3 Scaling properties of dust emission 

Some considerations must be made in order to enable the derivation of various physical quan- 
tities for each source from the comparison of the model expectation with the observational 
data. 

All models are computed assuming a luminosity of the primary power source of L^Giv _ 
Xq46 gj,g Both the intensity of dust emission and the torus sizes are strongly related to 
this value, as clearly indicated by Eq. about the relationship between the mininum radius 
and the AGN luminosity. The problem is that of finding how the solutions of the radiative 
transfer equation depend on the primary power source luminosity and on the size of the 
emitting region, and how we can scale the spectral solutions from the assumed L^'^'^ = 10^^ 
erg to the value require d by the target source. 



Rowan- RobinsonI (jl980|) notes that the radiative transfer problem has scaling properties 



that allow a reduction of the number of indepe ndent input parameters. These properties were 



first derived for spherical-shell geometries, but 



Ivezic fc Elitzurl ()l997l ) showed that they are 



a general property of dust emission in any arbitrary geometry. They showed in particular 
that main pivotal parameter in the dust emission modelling is the overall optical depth. So, 
as long as the optical depth remains the same, the size of the system can be scaled up or 
down arbitrarily, without affecting the radiative transfer solution. Obviously, the scale size 
of the torus distribution determines the value of the dust mass. 

Then, assuming a given value for the equatorial optical depth, if we change the source 
luminosity and size, we have the following scaling relations: 

Tocp-R^ pocR-^ ^ ^ocR-^ ^ M (xR^ (14) 

R-^ 

where R is the torus radius and M is the total dust mass. On the other hand, the torus 
luminosity is proportional to the grain number density and the grain emissivity e, that is 
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L oc pe/m oc M, where m is the average grain mass. Altogether, total dust mass and torus 
luminosity scale in the same way, like the square of the linear size. 

The mass of dust for a geometrical configuration of the torus is simply calculated as the 
sum of the mass of all the grains contained in each volume element, summed over all the 
elements of the grid. The mass of a dust grain is calculated from the typical density of the 
corresponding grain specie, by means of: 

4 

"^i = • af • p (15) 

where the i index is meant to represent grains with given size and p is the mass density 
of the specie, which was assumed to be 2.26 a nd 3.50 gr/cm^ for graphite and silicate 
respectively (see e.g. IWeingartner &: Draindl2001[ ). 

To summarize, the parameters of the torus model include the torus opening angle G, 
the ratio Rmax/Rmin of the outer-to-inner radii, the coefficients P and 7 that enter the 
density function, and the equatorial optical depth measured at 9.7 pm. Other parameters, 
like for example those describing the distribution function for the sizes of dust grains or their 
maximum and minimum values, were kept constant in order to avoid further complicating 
the problem. 

A final quantity to consider is the torus viewing angle \E': apart from specific cases, we 
will henceforth assume that it is either \E' = 0° for a polar view needed to represent type-1 
AGN, or \I' = 90° for an equatorial view needed for type-2 AGN. 



4 MODELING THE SED OF ACTIVE GALACTIC NUCLEI 



In the present Section we exploit our AGN IR model to study the observed broad band 
SEDs of a variety of mostly nearby objects, known to host active nuclei, selected from the 
NASA/IPAC Extragalactic Database (NED). We confine the study to sources with numer- 
ous (at least ten independent) available photometric datapoints in order to ensure a good 
sampling of their SEDs. 

A particularly interesting class of active galaxies are the Ultra-Luminous and Hyper- 
Luminous InfraRed Galaxies (ULIRGs and HyLIRGs, respectively). The studies of ULIRGs 
by IRAS and later by ISO often revealed the simultaneous presence of starburst and sig- 
nificant AGN dust components, whose r espective contribut i ons range from a s r nall fraction 
up to the totality of the IR SED (e.g 



2003; 



Prouton et al. 



Genzel et al 



199 



Verma et al 



2002; 



Berta et al 



2004^ . The IR emission of a starburst typically peaks at wavelengths 
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in th e range between ~ 50 and 100 \im (see e.g. 



Rowan-Robinson et al. 



1984; 



Soifer et al. 



1984^ . while the peak of AGN emission occurs at shorter wavelengths, somewhere between 



3 and 40 jim. Hence, given the fact that AGN and starburst phenomena often happen con- 
comitantly, and that in nearby AGN the host galaxy component can make an important 
contribution to the optical and near-IR spectrum, all three components have been included 
in our spectral modeling. For the cold dust component, the major contributor to the emis- 



sion at wave 



Farrah et al, 



ength s longer than ~ 50 pm, we chose to use starburst templates (see e.g. 
2OO3I ). More precisely, M82 represents a typical starburst IR emission while 



Arp 220 is adopted as representative of a very extinguished starburst. Additionally, NGC 
1482, NGC 4102, NGC 5253 and NGC 7714, are used in order to widen the choice of spectral 
shapes. The main differences between these templates are the intensity of the PAH features 
(between ~ 6 and ~ 15 jam), the depth of the silicate feature at 9.7 pm and the width, 
intensity and peak wavelength of the IR bump. 

A more exhaustive app roach would require to provide als o a physical description of the 
starburst component (e.g. lEfstathiou fc SiebenmorgenI 120051 ) . which however is far beyond 
the scope of this work. Our simplified approach of using instead observational templates of 
local starbursts is in line with the current work aimed at modelling the AGN dusty tori 
and deriving values for the best-fit parameters from comparison with the data: as explained 
below, starburst templates are added only when a torus model fails to provide an acceptable 
description of the observational SED. Concerning the starburst contribution, it is mainly 
determined by the far-IR datapoints, predominantly from IRAS 60 and 100 jim measure- 
ments, which are prone to considerable uncertainities. For this component we limit our 
analysis to just an inference of the bolometric IR luminosity. Finally, the stellar component, 
when needed, is modelled using opt ical model spectra of ULIRGs in the wavelength range 



from 0.3 to ~ 5 pm , as discussed bv iBerta et al. 



2OO3I and Fritz et al. (in preparation). Note 



that we left free the relative normalization of this optical/near-IR component and the far-IR 
spectrum of the host galaxy, given the extremely complex physical relation of the two. 



4.1 "Prototypical" sources 

In order to demonstrate the accuracy and reliability of our spectral modelling, we have 
studied in some detail the observed SEDs and IR spectra of three well known nearby objects. 
One of the three objects, Mrk 231, was chosen because it represents a type-1 AGN-dominated 
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Figure 10. The observed Spectral Energy Distribution of the Circinus galaxy (red triangles) including the 2.5 — 45 jam ISO SWS 
spectrum (thick cyan continuous line). The emission has been modelled using the NGC 1482 IR template spectrum (dashed 
line) plus a highly obscured AGN torus (thin line). The torus has an optical depth t(9.7) = 8.0, an aperture angle of 140° 
and a outer to inner radii ratio of 60, and an equatorial view. The bolometric AGN luminosity turns out to be ~ 1.78 X 10^* 
erg s~^. The density profile of dust follows a power-law radial profile with ind ex = —1 and ■y = 6.0. The integrated stellar 
emission component was taken from the optical spectrum of IRAS 19254-7245 iBerta et al.ll2o63 : Fritz et al. in prep.). 



ULIRG, while NGC 1068 was chosen as a prototype type-2 AGN-dominated source. The 
Circinus galaxy was finally considered as a composite object. All three are very well studied, 
given their proximity and brightness. 

Spectroscopic data in the mid-IR from ISO (SWS) were available for the three sources, 
imposing strong constraints to the torus characteristics. For this reason we built other ded- 
icated models with slightly different characteristics with respect to our reference grid (see 
Section IXT|) . to try to fit both the observed photometry and the mid-IR spectrum. 



4.1.1 The Circinus galaxy 



This very nearby (~ 4 Mpc) source displays characteristi cs typical o: 



like for example the very in tense forbidden o ptical lines (|01iva et al 



component of Ha emission (jOliva et al 



a Sey fert 2 object. 



1994^ and a broad 



19981 ). Moreover, indications of current star for- 



mation act ivity close to the ce ntral region of this spiral galaxy come from near-IR ob- 



servations (IMaio 



(Elmouttie et al 



ino et al. 



199 



19981 ) and from the extended morphology of the Ha emission 
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Figure 11. The SED of the Seyfert 2 galaxy NGC 1068, including the ISO combined LWS+SWS spectrum (thick magenta 
line) taken from Genzel & Cesarsky (2000). The model includes the contribution of the IR template of the starburst NGC 
7714. The assumed stellar spectrum is obtained as a composition of two models of the ULIRG IRAS 20551-4250 (Fritz et al. in 
preparation). The torus model has an aperture angle of 160°, an outer-to-inner radius ratio of 20 and a distribution of the dust 
density varying both in the radial and in the altitude coordinate (/3 = —1.0 and 7 = 6.0), with a value of the optical depth at 
9.7 pm of 8. The total luminosity is 3.9 X 10"*^ erg and the outer radius ~ 16.4 pc. 



The photometric datapoints availa jple for this source from the optical to the far-IR and 
the ISO SWS spectrum (2.5 to 45 pm; Moorwood et aLl[l996 . see Fig.fTUj) strongly constrain 
the AGN contribution. A large value for the equatorial optical depth [t(9.7) = 8.0, corre- 
sponding to Ay ~ 170 magnitudes] is required in order to depress dust emission between 2 
and 10 pm to match the observed spectrum, and a small values for the torus size (internal 
and external boundaries of 0.2 and 12 pc, with an aperture angle of 140°) is indicated to 
reproduce the longer wavelength part of the spectrum. The estimated mass of torus dust is 
~ 380 Mq, a low value ensuing our small best-fit torus size. 

Our torus model compone nt has a fla r ed dis c geometry and its main features are quite 



similar to those suggested by 



Ruiz et al. 



( 20011 ): small torus size (the parameter h/R^ax-, 
where h is the vertical size, is ~ 0.98 for this model) and a high visual extinction. The torus 
sizes of our solution are fairly consistent with their results. The most important difference 



concerns the relativ e staburst-AGN 
close to unity, while iMoorwood et al 



R contribution: 



T^iii7- et a1. 



( 20011 ) find this ratio to be 
( 19961 ) estimate the ~ 90 per cent of the total IR lumi- 



nosity to be ascribed to the AGN component. Our best-fit value for the AGN contribution 
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Figure 12. IR SED and spectrum of Mrk 231 versus model fit. To the far IR part of the SED the spectral template of NGC 
5253 has been added. This type-1 AGN has a torus with optical depth of 10, outer radius of 300 X Rmin with Rmin = 0.69 
pc and total torus aperture angle of 140° . The line-of-sight angle with respect to the equatorial plane of the torus is ^ = 25° , 
hence grazing the upper torus edge. 



is ~ 30 per cent, much closer to that proposed by iRowan-R.obinson fc CrawfordI ()1989[ ). Dif- 
ferences in the resuhs might be due to aperture correction effects, since the AGN emission 
i s concentrated in the inner galaxy bulge, while the starburst component is more extendend 



( Maiolino et al 



19981) . 



4.1.2 NGC 1068 



NGC 1068 is a very well studied prototype Seyfert 2 galaxy and the first one against whic h 
the Unified AGN Model was tested (see e.g. lAntonuccilll984l and lAntonucci fc Milled 1 198 51 ). 



Among oth ers, studies o f the p utative dusty toru s in this source included the IR polarized 



imaging by 



Young et al. 



Efs 



athiou et al. 



( 19951 ) ■ and recent mid- 



( 199fi[ ). the modeling b y 
IR interferometric observations by (jJaffe et al.l |200^ spatially resolving the dusty torus 
surrounding the active nucleus. 

The torus model that better fits the data (see Fig. ITT|l is one with an outer-to-inner 
radial ratio of 20 which, rescaled to fit the observed datapoints, gives a value of about 17.4 
pc for the external radius for a luminosity of ~ 3.9 x 10^^ erg s~^ of the central source. 
An optical depth of t(9.7) = 8.0 along the equatorial plane, together with a non-constant 
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density profile (/3 and 7 parameters for tlie densty law are —1.0 and 6.0 respectively), yield 
a silicate feature which is only slightly in absorption, as seen both in the observed SED and 
in the ISO spectrum. The derived optical extinction for such a configuration of the torus 
is Ay ~ 170, corresponding to a hydrogen column density of 6.98 x 10^^ cm~^ and a total 
mass of dust of 1.26 x 10^ Mq. 

The colder component of dust emission is well reproduced by adding the contribution of 
the IR template of the starburst galaxy NGC 1482, contributing some 70 per cent of the total 
IR luminosity. Our results on the AGN torus characteristics are in overall good agreement 
with previ ous studies , but i ndicate a larger torus than that inferred by the interferometric 
studies of jJaffe et al.l (|2004^ . who reveals the presence of dust at 320 K in a structure of 
3.4 pc in diameter. Anyway, in the model adopted to reproduce NGC 1068 emission in the 
mid-infrared, the dust reaches such a temperature within a radius of ~ 8 pc. Considering 
that the density law in this particular case takes the form: 

p(r,^) =a-r-^ •e-^^l^"^^^)! (16) 

we can easily see that the dust density dro ps in such a way t hat at 8 pc it is ~ 5 times lower 



than that at ~ 1.7 pc (radius reported by 
marginal agreement with the obsevations. 



■Taffe et al 



20041 ). Our findings are, therefore, in 



4.1.3 Mrk 231 



This object is one of the most luminous ULIRGs in the local universe fISoifer et al 



198J, 



and also one of the most studied at all wavelength. Both ROSAT and A SCA observations 



found clear evidenc e of a starburst activity combined to that of the AGN (jimanishi fc Ueno 



199fl : 



Turner 



Braito et al. 



199!^ ). The analysis based on Beppo-SAX and XMM-Newton performed by 



( 20041 1 indicates the presence of a very thick screen with an hydrogen column 
density of ~ 2 x 10^^ cm~^ covering the primary X-ray emission. 

The torus model is very strictly constrained by the observational data in the 1 to 5 pm 
interval. To reproduce the observed optical-to-IR SED (shown in Fig. IT^ we need a direct 
view of the central power source. The observed SED requires an additional contribution by 
colder dust at A > 10 pm, which we reproduce using the NGC 5253 IR SED. 

The torus model has t(9.7) ~ 10.0 (corresponding to Ay — 210) and a mass of dust con- 
fined in the molecular torus of about 2.7 x 10^ Mq, Rmax/Rmin = 300 and the torus aperture 
angle G = 140°. Such a high value for the optical depth implies a hydrogen column density 
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of ~ 9 X 10^^ cm~^, along the equatorial plane of the torus. The model AGN bolometric 
luminosity is 2.8 x 10^^ erg s~^, corresponding to an outer radius of ~ 200 pc. 

Our fit (whose parameters are summarized in Table 2) implies an AGN contribution to 
the IR flux of about 27 per cent, two thirds of the emiss i on be ing attributed to the host 



galaxy. This is consistent with the results by 



Braito et al. 



( 20041 ) ■ that found the starburst 



luminosity to be significantly higher than that of the AGN from a comparison of the X-ray 
and the bolometric flux and assuming a standard type-1 AGN emission. 



4.2 Further Applications 

In this Section we study the observed broad band SEDs of a sample of 40 type-1 Seyfert 
galaxies and quasars and 15 type-2 Seyferts, following the same procedure as in the previous 
Section but making use of our standard grid of models (see Section IXTj) . The aim is to obtain 
information about physical properties of the circum-nuclear dust distributions in a sample 
of local AGNs, as well as on the relative contribution of the host galaxy. 

In this exercise, we have taken particular care in selecting photometic measurements 
within homogeneous apertures, as much as possible. An example is given by our fit of the 
type-1 AGN NGC 7469 (Fig. IT^ . for which photometric data with an aperture including 
only the very central regions are used to constrain the intensity of the AGN component, 
while the total intensity fluxes include also the host galaxy contribution. 



4.2.1 Type-1 AGN 

A total of 40 type-1 AGN and quasars, with a good coverage of photometric datapoints, have 
been fitted with a combination of stellar, IR starburst and AGN emissions. The objects were 
chosen from various AGN samples with the only requirement to have a sufficient number of 
observational datapoints. Our first attempt was to reproduce the whole optical-IR observed 
SEDs with a pure AGN emission, trying all models in the grid and assessing their merit 
with the test. For those objects with unacceptably large obvious far-IR excess we 

added starburst and normal galaxy components in the far-IR and optical, respectively. 

In all cases the combination of the three emission components reproduced well all the 
observed SEDs. Figs. ^1 and El show the observed (red triangles) and model (black solid 
line) SEDs. The figures also report the various components of the model: the IR starburst 



A Dust Emission Model of AGN 27 



0.1 1 10 100 0.1 



10 100 0.1 1 10 100 



1 

0.1 
10-2 
10-^ 

1 

0.1 

10-2 

10-^ 

1 

0.1 

10-2 

^ 10-^ 

Ul 10-2 
10-^ 

1 

0.1 

10-2 

10-^ 

1 

0.1 

10-2 

10-^ 

1 

0.1 

10-2 

10-^ 




= 1 1 1 1 Mll| 1 1 1 1 1 lllj 1 1 1 1 1 Mlj 1 1 1 M MJ 

-IZwl _,^**x 

\ \^ \ 1 


= 1 1 1 1 M M| 1 1 1 M Ml| 1 1 1 1 Mllj 1 1 1 1 Mil 

:|IZw1 

r 1 1 1 Mill — 1 1 ii/iiif 1 1 mm — i i mm 


: 1 1 1 1 Mllj 1 1 1 1 M M| 1 1 1 1 M M| 1 1 1 1 1 1 Id 

:VIIZw244 : 

r 1 riTmif'^f^^ 


r 1 1 1 iiiii| — 1 1 mill — 1 1 mill — i ii iiin 
:Mrk009 = 
- ¥^ \ - 

" r*y^ \ " 


:Mrk010 : 
■ f^— \ 

r 1 1 1 Mill — 1 1 iiml|l 1 1 mm — i i iiira 


:Mrk079 : 
1 \ ! 

' V 

r ^ ^ \, \ 

\ : 


r 1 1 1 iiiii| — 1 1 mill — 1 1 mill — i ii lira 
^Mrk279 - '■- 

- jJr ^ 

r 1 1 1 iiiii| — 1 1 iifiif 1 1 mill — 1 II lira 


^Mrk290 

• 

I '-\\ 


r 1 1 mill — 1 1 1 iiiii| — 1 1 1 iiiii| — 1 1 mm 
^Mrk304 '-- 

r^'i iiiiii| — 1 1 1 iiiii| — 1 1 1 iiiii| — 1 1 mm 


^Mrk335 = 

■ ^' V 

r 1 1 1 iiiii| — 1 1 iii/i'i 1 '1 1 mill — 1 II lira 


?i 1 1 Mill — 1 1 1 iiiii| — 1 1 mm — 1 1 1 1 IB 
:Mrk359 

: r \ '■ 
r 1 1 1 Mill — 1 1 iiiiif 1 1 mm — i i iiira 


:Mrk376 . : 
/Y / \ \ n 

^ ^ N 


:Mrk463 

r /</ \ . 

y ^ \ \: 


:Mrk478 

r 1 1 1 Mill — 1 1 iiiyii|' 1 1 mm — i i mm 


rMi mill — 1 1 1 iiiii| — 1 1 1 iiiii| — 1 1 mm 
:Mrk486 

i ^ 


r 1 1 1 iiiii| — 1 1 mill — 1 1 mm — i 1 1 lira 
^Mrk506 : 

^ / \ ■ 

>^ 1 M 1 X 


:Mrk509 : 
^ \ 1 


- I I mill — 1 1 1 iiiii| — 1 1 1 iiiii| 1 1 mm 
:Mrk704 : 

r ^^^^ 1 

/ \ ! 

- ^ / \ 

y \ : 


; 1 1 1 lllil| — 1 1 1 lllll| — 1 1 llllll| — 1 1 1 lira 

:Mrk734 : 

r ..,^^^^\^ 

— 1 .1-1 1 ml 1 1 1 1 1 1 1 1 1 m 


rill iiiii| — 1 1 1 iiiii| — 1 1 iiiiii| — 1 1 mm 
:Mrk771 

- \^ 
1 1 1 1 1 ' I'll"* 


r 1 1 iiiiii| — 1 1 1 iiiii| — 1 1 1 iiiii| — 1 1 1 Ilk 
:Mrk876 ^ 



10-2 

10-^ 

1 

0.1 

10-2 

10-^ 

1 

0.1 

10-2 

10-^ 

1 

0.1 

10-2 

10-^ 

1 

0.1 

10-2 

10-^ 

1 

0.1 

10-2 

10-^ 

1 

0.1 

10-2 

10-^ 



1 1 10 100 0.1 1 10 100 0.1 1 10 100 



Figure 13. Comparison of our best-fit models with observational data for the type-1 AGN. Red triangles are the observed data 
points; the solid black line is the total model SED; green line is the AGN emission; the blue dashed line denotes the star burst 
emission; stellar emission by the host galaxy is shown as dashed lines. 
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Figure 14. Comparison of our best-fit models with observational data for the type-1 AGN. Meaning of symbols and line is 
the same as for Fig. 1131 



A Dust Emission Model of AGN 29 



Obj 


z 


LAGiV 




^Max 


e 


/3 


7 


r(9.7) 


Av 






AGN 








[lO^^erg s-l] 


[pc] 


[pc] 


O 








[mag] 


cm' 


-2 


Mq 


[%] 


O 


3C249.1 


0.31120 


0.89 


1.23 


36.8 


100 


-1.0 


6.0 


10.0 


214.5 


8.68 


X 


Iq23 


6.85 X 10^ 


72.2 





IRAS 13349 


0.10764 


0.79 


1.16 


116 


140 








1.0 


21.6 


8.79 


X 


1022 


1.13 X 10^ 


100 





IRAS 13517 


0.08820 


0.32 


1.29 


129 


100 


-1.0 


6.0 


10.0 


214.5 


8.68 


X 


Iq23 


6.20 X 10'' 


29.6 


20 


IRAS 16347 


1.33400 


7.94 


3.66 


1100 


140 








2.0 


42.7 


1.73 


X 


Iq23 


1.99 X 10^ 


100 


20 


[HB89] 1543+489 


0.40000 


0.56 


0.98 


29.2 


140 


-1 





10.0 


214.5 


8.68 


X 


1023 


3.12 X lO'' 


55.3 





[HB89] 1700+518 


0.29200 


0.79 


1.69 


50.7 


140 


-1 





6.0 


129.1 


5.23 


X 




5.42 X 10'' 


68.1 





IZwl 


0.06114 


0.25 


0.65 


65.2 


100 








2.0 


43.5 


1.76 


X 


1023 


5.80 X lO"* 


35.0 





IIZwl 


0.05434 


0.14 


0.49 


147 


100 


-0.5 


6.0 


1.0 


21.9 


8.94 


X 


1022 


1.47 X 10^ 


14.3 





VIIZw244 


0.13100 


0.14 


0.49 


14.7 


140 








6.0 


129.0 


5.23 


X 


1q23 


1.16 X lO" 


76.2 





Mrk 009 


0.03987 


0.07 


0.35 


104 


100 


-0.5 





1.0 


21.9 


8.94 


X 


1022 


4.16 X lO"* 


47.6 





Mrk 010 


0.02925 


0.03 


0.21 


61.8 


100 


-0.5 





3.0 


64.9 


2.64 


X 


Iq23 


4.36 X 10" 


53.3 





Mrk 079 


0.02219 


0.01 


0.15 


4.5 


140 


-1 





10.0 


214.5 


8.68 


X 


Iq23 


7.30 X 102 


46.5 





Mrk 231 


0.04217 


0.28 


0.69 


207 


140 


-0.5 





6.0 


128.8 


5.22 


X 


1023 


1.34 X 10^ 


27.5 





Mrk 279 


0.03045 


0.03 


0.23 


69.4 


100 


-0.5 





3.0 


64.9 


2.64 


X 


1023 


5.48 X 10'' 


49.1 





Mrk 290 


0.02958 


0.01 


0.15 


4.4 


140 


-0.5 





6.0 


129.43 


5.25 


X 


1023 


7.01 X 102 


53.8 





Mrk 304 


0.06576 


0.15 


0.51 


15.2 


100 


-1.0 


6.0 


10.0 


214.5 


8.67 


X 


Iq23 


1.16 X 10^ 


64.5 


20 


Mrk 335 


0.02578 


0.18 


0.55 


16.4 


100 


-0.5 


6.0 


6.0 


129.5 


5.26 


X 


Iq23 


1.38 X 10^ 


88.3 





Mrk 359 


0.01738 


0.01 


0.10 


31.0 


140 








6.0 


127.9 


5.16 


X 


Iq23 


4.69 X lO"* 


59.8 





Mrk 376 


0.05598 


0.06 


0.33 


9.8 


140 


-1 





10.0 


214.5 


8.68 


X 


Iq23 


3.49 X 10^ 


57.6 





Mrk 463 


0.05035 


0.14 


0.49 


148 


140 








6.0 


127.9 


5.16 


X 


Iq23 


1.08 X 10'^ 


71.0 





Mrk 478 


0.07905 


0.32 


0.73 


21.9 


100 


-1 


6.0 


10.0 


214.5 


8.68 


X 


1q23 


2.43 X 10^ 


50.8 





Mrk 486 


0.03893 


0.08 


0.37 


11.0 


140 


-0.5 


6.0 


10.0 


214.3 


8.67 


X 


Iq23 


1.03 X 10^ 


100 





Mrk 506 


0.04303 


0.04 


0.26 


25.9 


100 


0.0 





3.0 


64.5 


2.62 


X 


1023 


1.18 X 10'' 


86.4 





Mrk 509 


0.03440 


0.09 


0.40 


11.9 


140 


-0.5 





6.0 


129.4 


5.25 


X 


Iq23 


5.20 X 10^ 


65.5 





Mrk 704 


0.02923 


0.04 


0.25 


7.4 


140 


-1 





10.0 


214.5 


8.68 


X 


1023 


1.97 X 10^ 


90.0 





Mrk 734 


0.02923 


0.03 


0.21 


61.8 


140 








2.0 


42.7 


1.73 


X 


1023 


6.29 X 10'' 


51.5 





Mrk 771 


0.06301 


0.02 


0.20 


5.8 


140 


-0.5 





10.0 


214.3 


8.67 


X 


Iq23 


6.51 X 10^ 


41.9 





Mrk 876 


0.12900 


0.28 


1.22 


122 


100 


-1 


6.0 


10.0 


213.7 


8.64 


X 


1023 


5.52 X 10'' 


23.5 





Mrk 1014 


0.16301 


0.89 


1.23 


123 


140 


-1 


6.0 


10.0 


213.6 


8.63 


X 


Iq23 


6.51 X 10" 


33.2 





Mrk 1383 


0.08647 


0.26 


0.67 


20.0 


140 


-0.5 





10.0 


214.3 


8.67 


X 


1023 


2.42 X 10" 


79.7 





Mrk 1513 


0.06298 


0.28 


1.21 


121 


100 


-1 


6.0 


10.0 


213.6 


8.63 


X 


1023 


5.40 X 10" 


70.7 





NGC 3516 


0.00884 


0.01 


0.12 


3.7 


100 


-1 


6.0 


10.0 


214.5 


8.68 


X 


1023 


7.0 X 10^ 


32.8 





NGC 5548 


0.01717 


0.03 


0.23 


6.6 


140 


-0.5 





10.0 


214.3 


8.67 


X 


1023 


2.60 X 10^ 


67.5 





NGC 7469 


0.01632 


0.03 


0.23 


6.9 


100 


-1 





10.0 


214.5 


8.68 


X 


1023 


1.25 X 10^ 


8.8 





PG0804+761 


0.09990 


0.85 


1.20 


36.0 


100 


-1.0 


6.0 


10.0 


214.5 


8.68 


X 


1023 


6.54 X 10=* 


92.1 





PG0906+484 


0.11800 


0.16 


0.52 


15.5 


100 


0.0 





10.0 


214.0 


8.65 


X 


Iq23 


6.27 X 10^ 


28.7 





PGlOOl+054 


0.16100 


0.96 


1.27 


38.1 


100 


-0.5 


6.0 


6.0 


214.3 


8.67 


X 


1023 


4.45 X 103 


72.5 





PG1211+143 


0.08090 


0.59 


0.77 


23.0 


100 


-1 


6.0 


10.0 


214.5 


8.68 


X 


Iq23 


4.53 X 10^ 


63.8 





PG1411+442 


0.08960 


0.28 


0.69 


20.7 


140 


-0.5 





10.0 


214.3 


8.67 


X 


1q23 


2.59 X 10" 


100 





PG2112+059 


0.46600 


2.51 


2.06 


206 


100 


-0.5 





10.0 


214.2 


8.67 


X 


1023 


1.66 X 10'^ 


100 





UGC00545 


0.06114 


0.43 


1.50 


150 


100 


-1 


6.0 


10.0 


213.7 


8.64 


X 


Iq23 


8.36 X 10" 


29.5 






Table 2. Best-fit physical parameters of the torus model for all type-1 AGN. In the order, object name, redshift, AGN 
luminosity, inner and outer boundary radius, the full-opening angle of the torus, the radial and angular density dependency, 
the optical depth at 9.7 jam, the extinction Av, the hydrogen column density J^h, the mass of dust and the percentage 
contribution of the AGN. The latter is computed in the range 5 + 1000 jam. The viewing angle 't, measured with respect to the 
2 axis, is in general taken to be 0°, since the models show in general very little dependence on this parameter for type-1 AGN. 



and AGN contributions on blue dashed and green solid lines, respectively, while the stellar 
emission is shown as dotted lines. 

Table 121 summarizes the physical parameters derived for the 41 type-1 AGN under inves- 
tigation (including Mrk 231). In the order, the columns report the object name, redshift z, 
the luminosity of the central power source, the inner and outer torus radii, the full-opening 
angle of the torus, the radial and angular dust density dependences (parameters (3 and 7, 
respectively), the optical depth at 9.7 jam, the extinction Ay, the hydrogen column den- 
sity Nh, the mass of dust and the percentage contribution of the AGN to the IR bolometric 
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Table 3. Summary of the main features of the torus model used to reproduce the SED of the 17 type-2 objects (including 
Circinus and NGC 1068). /3 and 7 are the parameters entering the density law, is the full-opening angle of the torus, t(9.7), 
Nh and Ay refer to the equatorial plane of the torus. 



emission. Note that t(9.7), Ay and Nh are measured between the inner and the outer radius 
in the torus equatorial plane. 

For the majority of the objects (27 out of 40), the contribution of the AGN torus to the 
total IR emission is larger than 50 per cent. The SEDs of five objects in particular, Mrk 
486, IRAS13349, IRAS16347, PG1411+442 and PG2112+059 can be entirely reproduced 
using a pure AGN model, fitting of their observed SED does not require significant cold 
dust emission. 

As for the short-wavelength optical-UV part of the spectrum, there is a remarkable 
complementarity of the spectral shapes for the AGN and host galaxy: the local minima of one 
component correspond to the maxima of the other, including the 9.7 pm emission/ absorption. 
This coincidental fact allows us to assess with some precision the contribution by the host 
galaxies to optical spectrum. Quite often in these moderate luminosity Seyfert-1 galaxies 
the latter appears to be important (20 cases out of 40). 

Our model solutions for the type-1 objects in this sample tends to favour dust tori in 
which the density has a radial gradient: 32 of the 40 sample sources require a power-law 
density profile decreasing with radius (i.e. /3 = —0.5 or P — —1.0), and 6 of them also favour 
an angular dependence (7 = 6.0). 

Typical torus sizes are Rmax < 100 pc. Only for 9 objects an outer radius barely above 
100 pc is indicated, with the exception of the highest rcdshift one, IRAS 16347+7037, that 
has a very luminous central source (^^ 8 x 10"^^ erg s~^) and whose external radius is found 
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Figure 



15. Spectral models for the type-2 AGN. Colour and line coding is the same as in the previous figures. 
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to be around 1 kpc. As for the torus geometry, it is remarkable that the majority of our 
sample sources are well fitted by an high total aperture angle (0 = 140°), implying an high 
dust covering factor. 

For the large majority of cases, the values of the mass of dust range from few tens to 
few tens of thousands solar masses, orders of magnitudes lower than the typical masses of 
supermassive black hole. Striking exception is again the mass of dust computed for IRAS 
16347+7037 found to reach ~ 10^ Mq. 

Most important, the best-fit values of the optical depth at 9.7 \im were found to be larger 
than 1 for almost all the objects, corresponding to visual extinction values larger than 22 
magnitudes and hydrogen column densities in the torus that must exceed ~ 10^^ cm~^. 

4.2.2 Type- 2 AGN 

The same fitting procedure has been used to analyse the observed broad band SEDs of a 
small sample of 15 type-2 AGN. Table 121 summarizes the main best-fit parameters (including 
those for Circinus and NGC 1068). The results of our SED modelling (shown in Fig. ITH|l 
differ in some aspects from those for type-1 objects. 

Again, all observed SEDs can be very well reproduced by combining stellar, starburst 
and AGN emission. Note that, in this case, the optical and near-IR emission is entirely due 
to the host galaxy's stellar contribution, since the AGN emission is totally suppressed at 
wavelengths shorter than 2 — 3 pm. None of the type-2 objects can be reproduced with torus 
emission alone: they all require a substantial starburst component. An Arp 220 starburst 
template better reproduced the far-IR data in eight out of 15 cases. 

As shown in the left panel of Fig.^J summarizing a comparison of our best-fit parameters 
for the two AGN populations, the average AGN contribution to the bolometric IR luminosity 
is lower than that estimated for the type-1 objects. We should caution, however, that the two 
subsamples are distributed quite differently in distance: the redshifts for the type-1 objects 
(appearing in the same Fig. fT6|) are typically higher than for the type-2, due to their easier 
spectroscopic identification. This bias obviously affects the bolometric AGN luminosity and 
dust torus mass (Fig. IT7|l . while it is unclear at this stage how much it might influence the 
estimate of the AGN fractional contribution. 

It is remarkable that, in spite of this bias, the distributions of the main physical and 
geometrical parameters for the two classes of AGN are fairly consistent, in particular the 
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Figure 16. Left panel: fractional contribution of the AGN component to the IR bolometric (5 — lOOOym) luminosity versus the 
redshift, for type-1 AGN (red 4-leg stars) and type-2 (green triangles) AGN. Right panel: summary of the main characteristics 
of the torus best- fit parametric solutions, with a comparison of type-1 (red continuous line) and type-2 (green dashed line) 
sources. 
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Figure 17. The distribution of the values of dust masses for type-1 and type-2 sources does not seem to highlight deep 
differences between the two classes, while distribution of the IR luminosity that is ascribed to the AGN component is more 
peaked towards the low luminosity part of the plot. This effect is due to the fact that the star burst component is more important 
in type-2 sources. 



optical depth, with a median value of t(9.7) of 6.0 for the type-2, while it is 10.0 for the type- 
1 sample. In any case, the results should be seen with caution due to the in- homogeneity, 
incompleteness and random character of our sample. 
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5 DISCUSSION AND CONCLUSIONS 

This work describes an improved model for the emission of dusty tori around AGN, including 
features already occasionally discussed in the literature and paying particular attention to an 
accurate computation of the radiative transfer solution. The chosen torus geometry is simple 
but realistic, a flared torus defined by i ts inner and outer ra dii and the total opening angle. 



We adopt the dust grain distribution of iMathis et al.l ()l977f ). with the density depending on 



the torus radial coordinate and polar angle. The optical depth is computed in great detail, 
taking into account the different sublimation temperatures for silicate and graphite grains. 

The fact that grains with different sizes can have different sublimation temperatures may 
imply a significant complication, especially in the inner parts of the torus where the tem- 
perature of some dust grains can reach and exceed the sublimation limit. While this model 
accounts for the differences between graphite and silicate gr ains, we assume the sublimation 



tempe rature to be independent of the grain size. This follows 



Efstathiou fc Rowan-Robinson 



( 19941 ). who have shown that this effect is not of much importance especially for high den- 
sities and optical depths. 

Our model is able to reproduce with remarkable accuracy the observed broad-band SEDs 
of a variety of AGN, for which large numbers of photometric measurements were available 
in the literature. For three well-known nearby active galaxies (Circinus, NGC 1068 and Mrk 
231), the ISO mid-IR spectrum was also available, in which cases both the torus model and 
the starburst contributions were strongly constrained. 



5.1 Spectral fitting of type-1 AGNs 

A starburst component is usually required in order to fit the observed SED longward of 50 
pm, as in most cases a pure AGN component is unable to reproduce the observed width 
of the IR bump. From the objects of our sample, only four type-1 sources can be fitted by 
a pure AGN emission in the far-IR. The torus emission is in general predominant for such 
objects, while the starburst prevails in type-2 (see Fig. UHl left panel). 

In this work we pay particular attention to the behaviour of the 9.7 pm silicate feature, 
which in many previous publications was found to challenge the Unified Model (especially 
for type-1 objects). The problem was the prominence of the feature in absorption for type-2 
AGN and its absence in type-1 objects, as viewed pole-on. This feature arises from dust sit- 
uated in the innermost regions of the torus, as seen in Fig. El were the temperature is close 
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Figure 18. Model fit to the observed SED and the IR Spitzer spectrum (observed by Siebenmorgen et al. 2005) between 
5 — SOpm of the luminous QSO 3C249.1. The 9.7 jam silicate feature is clearly seen in emission with an intensity which very 
close to that of our model. The discrepancy in the wavelength of the observed peak is likely due to a different chemical 
composition of the silicate grain I Siebenmorgen et al. 2005; Hao et al. 2005). A non-constant density profile with 7 = 6 and 
/3 = —1.0 and a torus aperture angle of 100° arc used, while the external radius is ~ 35 pc. 

to the sublimation limit. Moving far away from the center the feature becomes gradually 
fainter untill it eventually turns into absorption if the optical depth is high enough. Our cur- 
rent results show that a combination of high optical depth, moderate torus radius, an either 
constant density profile or one radially increasing toward the center, and a proper consid- 
eration of the lower sublimation temperature of silicate grains, reduces the silicate emission 
or even turns it into slight absorption. Under these assumptions, a standard composition of 
the mixture of dust turns out to be entirely compatible with low or no silicate emission in 
our model. The addition of a significant starburst contribution with the silicate feature in 
absorption tends to further decrease the prominence of the 9.7 jim feature in type-1 AGN. 

Recent observations with the Spitzer IR space telescope have revealed , for the first time, 
evidence for the silicat e feature in moderate emission in a few type-1 AGN (jSiebenmorgen et al, 



Hao et al, 



20051 ) ■ with values of the 6*9.7 parameter close to unity. Indeed, the latter 
can be achieved with several combinations of our torus parameters (see Fig. E)). We present 
in Figs. ^1 and El a detailed comparison of IR Spitzer spectra, together with SED data in 
the optical and far-IR, of two quasars-1 (3C249.1 and PG0804+761) with evidence for the 
9.7 pm feature in emission. We see that our model provides an excellent fit of the spectral 
shape and broad-band data for the latter source. As for the former, the spectrum is also well 
reproduc ed, except for a slight shif t in wavelength probably due to a locally different grain 
mixture ( Siebenmorgen et al.ll2005 ). In conclusion, we find no need for particular geometries, 
like a tapered disc, or for modified distributions of the silicate grains (clumpy tori) in order 
to suppress the silicate feature in type-1 sources. 
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Figure 19. Comparison of the Spitzer spectrum iHao et al.ll20nfi^ and other SED datapoints with our model expectation for 
the typo-1 AGN PG0804+0761. The model is similar to that fitting 3C249.1 (see Fig.EHand parameters in Table 2). In this 
case our fit to the silicate peak is very good. 
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Figure 20. Comparison of predicted spectra for a torus aperture angle of 140° (continuous line), 100° (dotted line) and 60° 
(dashed line), for a model with the same characteristics as the one used in Fig. 1191 Datapoints as discussed in the caption to 

Fig. Eg 



5.2 Torus geometry and the Unified Model 

We now discuss results related to the Unified Scheme. Fig. (right panel) illustrates the 
best-fit values of the parameters reproducing the observed SED of type-1 and type-2 AGN. 
Non-constant dust density profiles, in particular a density decreasing with the distance from 
the center, are clearly favored. Furthermore, high values for the equatorial optical depths are 
required in almost all the sources, particularly to suppress the 9.7 jim emission. Interestingly, 
in agreement with the predictions of the Unified Scheme and despite the different average 
distances for the type-1 and type-2 samples, the distributions of the geometrical and physical 
torus parameters do not show significant variation between the two populations. 

Also the distribution of values for the mass of dust in the torus and the luminosity of 
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the primary power source presented in Fig. El do not show dramatic differences between 
type-1 and type-2. The shght tendency of type-1 objects to host more luminous AGN and 
more massive tori may be explained with the Malmquist bias introduced by the different 
average distances. More problematic is, instead, the interpretation of the higher fraction of 
the AGN contribution in type-1 sources, seen in Fig. (left panel). It is unclear whether 
this is still due to a distance-induced bias, or if it is reflecting an intrinsic differentiation of 
the two populations. 

Obviously, not only the silicate feature, but more in general the observed mid-IR spectral 
shape can effectively constrain the intrinsic properties of the circum-nuclear dust distribu- 
tion. An illustration of this is given in Fig. 1201 showing how the spectrum changes by reducing 
the torus full-opening angle 6 from 140° to 100° and 60°: the reduced amount of cold dust 
close to Rmax foT the 60° case makes the spectrum very fiat in flux density units from ~ 3 to 
~ 50 pm, whereas a much steeper one is produced by our more standard value of O = 140°. 

A further important implication of our work is that, even accounting for important 
starburst contributions to the longer- wavelength far-IR part, our code was essentially unable 
to obtain good spectral fits to both AGN classes by assuming relatively fiat tori with low 
values of the opening angle = 60°. Instead roughly half of the type-1 AGN and 70 per cent 
of the type-2 require almost cylindrical tori with O = 140°. For these, the corresponding 
covering factor is ~ 90 per cent, while the remaining sources are better fit with B = 100° 
models, corresponding to ~ 75 percent covering factor. 

Altogether the average covering factor for our local sample is close to 80 — 85 percent: 4 
to 5 obscured AGNs are then to be expected on average for each type-1 unobscured object 



[we predict very thick circum-nuclear material 



samples - including those selected in t 



in radio ( Lawrence 



1991 



le IR (e.g . 



but see also 



Maiolino 



Ay > 100 



Rus 



1 et al. 



. Since local and high-z A GN 



1993 



La Franca et al. 



20041 ) or 



19951 ) - do not typically reveal such a high 



incidence of obscured objects, it remains that a large number of completely obscured AGN 
are to be expected, undistinguishable from normal galaxies if observed in the optical/ near- 



IR or soft X-rays. This appears consistent with the results of 



Francesch ni et al 



a deep combined Spitzer/ Chandr a survey. As discussed there (see also 



Gandhi et al, 



mi 



Fabian 



Maiolino et al 



( 2005) from 



2m 



20031 ). an important fraction of the obscured quasars and AGN 



may escape identification, even with a good coverage of the mid- and far-IR spectrum, as 
given by the Spitzer data, or deep hard X-ray data. 
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5.3 Model degeneracies 

Inevitably, our solutions for the circum-nuclear dust structures in AGNs suffer some degree 
of degeneracy in the parameter space. In our case, the largest source of uncertainty is the 
limited spectral coverage, especially in the mid-IR. As already shown for the three proto- 
type objects analysed in Sects. 4.1.1, 4.1.2 and 4.1.3, the mid-IR spectrum turns out to be 
foundamental constraining the model parameters of AGN dust tori. The shape of the mid-IR 
(1 — 12 pm) continuum is ruled by the optical depth and the characteristics of the density 
law, while the shape of the far-IR SED (20 — 400 jim) is mainly determined by the torus 
size and the amount of cold dust. 

In order to study the degree of degeneracy of our fits, a detailed case-study of the 
following three objects was performed: PG1411+442, as representative of type-1 objects with 
no evident starburst contribution, Mrk 1513 as a type-1 source with significant contribution 
from starburst emission in the IR, and NGC 1614 as a representative of the type-2 sources. 

We have compared the observed SEDs for the three sources with all model spectra 
corresponding to the whole parameter set in Table ^ and then considered the solutions 
with a reduced below xlest + 3- The accretion luminosity turned out to be the best 
constrained quantity, especially in the type-1 sources where the primary continuum emission 
is observable. The scatter was found to be somewhat larger for the type-2 NGC 1614, in 
which case it varied between ~ 0.75 x 10"^^ and ~ 10^^ erg/s, the higher values corresponding 
to model solutions with smaller torus opening angles. 

The values of the outer radius were fairly constant for the best fit models in NGC 1614, 
but differed by up to a factor of 3 in the case of PG141H-442 (both Rmax/Rmin = 30 and 
100 gave acceptable fits). 

Similar results were obtained for the mass of dust. The most extreme case was again 
PG141 1+442, with a variation by up to a factor of 4, as a consequence of the Rmax/Rmin 
variation. On the other hand, the parameters of the dust density law (equatorial optical 
depth and the coefficients (3 and 7 of the density law) were found to only slightly vary 
among the best fits. The starburst template that better fitted the objects was Arp 220 and 
was kept fixed as no acceptable fits could be obtained with other templates. 

In conclusion, the available data are not so detailed to allow resolution of all degeneracies 
in the torus parameters, which we may understand if we consider the complexity of the 
objects under scrutiny and the fact that the torus and starburst emissions overlap smoothly 
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in the far-IR. However, our analysis has shown that significant constraints can be achieved 
about the characteristics of the circum-nuclear dust distributions in AGNs. 
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